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Numerical Determination of Characteristic Numbers J 

By W. E. Milne 

This paper represents a contribution to the problem of characteristic values and charac- 
teristic solutions of ordinary linear differential equations. The problem is conceived as a 
vibration problem in x and t. The partial differential equation is then approximated by a 
difference equation in both entries. The problem is now to find those frequencies in which 
a separation in x and t takes place. This is done by finding the roots of a trigonometric 
expansion of certain order. The method is applied to a number of interesting cases that 
illustrate various types of situations encountered in problems of physics and engineering. 



An improved method for obtaining the latent 
roots of a matrix has been devised by Cornelius 
Lanczos of this Bureau. He has shown further 
how this process may be used to secure the char- 
acteristic numbers (eigenwerte) belonging to a 
boundary value problem associated with an ordi- 
nary linear differential equation. The purpose of 
this note is to present a procedure for calculating 
characteristic numbers, based essentially on this 
method of Lanczos, but modified in such a way 
as to provide a simple numerical routine for the 
computation. A number of numerical examples, 
worked in full, illustrate the procedure and give 
an indication of the accuracy attained. The 
exposition is limited to the case of differential 
equations of the second order, but the method is 
capable of extension to cases of higher order. 

1. The problem. Let 



d 2 u 



du 



L(u)=P (x) %j+Pi(x) ;g+P«(*)tt. 

The problem before us is to find those characteris- 
tic values of X for which the differential system 



. du A , 

u-\-g , =0 at x=a 

u-\-G j- =0 at x=b 
ax 



(1) 



possesses nonzero solutions in the interval 
aSx^b. 

Problems of this type arise in many different 
ways in mathematical physics. For example, 
they may occur in connection with the heat 
equation 



dt 



=i<F), 



if we assume particular solutions of the form 

V=u(x)e-™, 

or in connection with the wave equation 
b 2 V 



1 The preparation of this paper was sponsored (in part) by the Office of Naval Research. 
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if we assume particular solutions of the form 

V=u(x) cos \t, 

or they may occur after separation of variables in 
more general partial differential equations. 

The important fact is that the characteristic 
numbers are completely determined by the sys- 
tem 1, regardless of the particular physical prob- 
lems from which 1 may have been derived. 
Consequently, in order to obtain numerical solu- 
tion we may assume a fictitious partial differential 
system which leads to (1). The most satisfactory 
fictitious system for our purpose appears to be 



b 2 V 



with boundary conditions 

bV 

V+g -§^=0 at x-- 

V+G^=0 atz= 



(2) 



-a, 
--b, 



(3) 



and initial conditions 

V=F(x) when i=0, 

=0 when t=0. 



bV 
bt = 



(4) 



Following standard procedure, we assume that 
the solution of (2, 3, and 4) is 



V=^c k u k (x) cos \ k t, 

k 



(5) 



where the X^ are the desired characteristic num- 
bers, the u k (x) are the characteristic functions 
(normalized in some definite manner) and the c k 
are constants depending on F(x) . 

2. Normalization of L(u). It will be assumed 
with respect to the coefficient functions in L(u) 
that for all values of x in the interval a^xSb 

(a) P (z)>0, 

(b) -T- P (x) exists and is continuous, 

(c) Pi(x) and P 2 (x) are continuous. 

We find it desirable to change the independent 
variable x to a new variable 5 by the relation 



dx 



VPoW 



=ds. 



This change of variable carries L(u) into what we 
shall call the " normal form" 



~'ds 2 



L(u)=^+p(s)^ + q(s)u } 



du 

} ds' 



where p(s) and q(s) are seen to be continuous in 
the interval a'tkstkb' corresponding to a^x^-b. 
It will be assumed from now on that L(u) has 
already been put in the normal form. 

3. The partial difference equation. Let the in- 
terval a^x^b be divided into n+1 equal sub- 
intervals of length 

t o—a 

71+1 

and let the value of V(x, t) at the point x=a-\-ih 
and at time t=jh, be designated by V^. Then 
in the differential equation 

we make the substitutions 



u 2 




h 2 


b 2 V 


v t+h 


,-2V tt +V t , 1 . i 


d 2 x 




h 2 


dV_ 


v i+1 . 


i-Vw 


dx 




2h 



and after simplification obtain the approximate 
difference equation 



hpi 



~~o~ H H-1.J" 



Vt-LA+VqiVv, (6) 



in which 



Pi=p(a+ih), q i= q(a+ih). 

If we assume the existence and continuity of the 
second derivatives of p(x) and q_(x), and suppose 
that the error of eq 6 is expanded in powers of h, 
we find that the leading term of the error is 



k 

12 



bV 



4 r b 2 V bV 

T(q"+pq' + q 2 )v], (7) 
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showing that the error of eq 6 is of fourth order 
in h. 

4. Derivation of the characteristic equation. For 
simplicity in the following explanation, let us 
suppose that in the boundary conditions g=G=0, 
so that V oj =0, F n+U =0. For £=0 (or j = 0) the 
values of V i0 are arbitrary (being equal to the 
arbitrary function F(x)). Hence we choose these 
in the most convenient manner, say Fio^O, 
V i0 =0, £>1. Then we have for V tj the array of 
values 



\ 





1 


2 


3 


. . . 


n 


w-fl 








v 10 
















1 





V n 


V n 













2 





Vu 


^22 


v 32 










3 





V\Z 


Vn 


^33 










n-1 





Kin— 1 


Vj»— 1 


V& n-1 




V n n- 1 





n 





V 2n 


V ln 


v 3n 




y nn 






(8) 



As just explained, the entries of the first row are 
arbitrary and are chosen as shown for convenience. 
To obtain the second row we use eq 6, together 
with the relation 

derived from the condition dV/dt=0 at t=0. All 
remaining rows are computed by eq 6. 

So far no restriction has been placed on the 
magnitude of the interval h. We now impose the 
condition that h is chosen small enough that 

I*Ap<I<i. 

for all values of i. 

It will then be seen from array 8, together with 
eq 6, that the entry in the principal diagonal in 
the second column is given by 

V 21 =m-ihp 2 )V 10 , 

whereas for subsequent columns the entries in the 
principal diagonal are 

V 32 =(l-hhp d )V 2U 

etc. 

Since by our choice Fio^O no entry in the 
principal diagonal is zero, and the ft-square matrix 
obtained from (8) by dropping the first and last 



columns and the last row is not singular. Then 
there will exist n constants A , A i} . . . , A n -i such 
that 

j^AjVij^O for*=l, 2, . . ., n, 

j = 

(where A n =l). 

Now from eq 5 we have 

n 

Vij=^Zc k u ki cos \ k hj, 

k=l 

where u ki =u k {a-\-ih) . Therefore, 

n n 

ScAiSi; cos \ k hj=0, i=l, 2, . . ., ft. 

fc-l ;=0 

It can be shown that these equations imply 

n 

y^.A< cos \ k hj=Q, k = l, 2, . . ., n. 

j = 

It is apparent, therefore, that the desired charac- 
teristic numbers \ k are given by 



X* = /i*/A, i=l, 2, 



ft, 



where ii Xi /x 2 , . . ., m« &re roots of the character- 
istic equation 



cos nn-\-A n _i cos (ft— 1)m+ ■ 



•-{-Aicos ii-{-A o =0. 
(9) 



The \ k thus found are correct for the difference 
equation with the given boundary conditions but 
in general are only approximations to the A* 
belonging to the differential system 1 . 

5. The characteristic functions. In order to 
obtain numerical values for the characteristic 
function u ki associated with \ k) we substitute 

Vij = U ki COS juj, 

in eq 6. Since 

V itj j rl +V i , j ^ 1 =2V ij cos n k , 

the result, after removal of the factor cos p k j, may 
be put in the form 

(l+%hpt) u k , w =(2 cos /xjt— h 2 g: t )u ki — 

(i— j*pO ^.i-i. (io) 

This is an ordinary difference equation of second 
order. We may choose any nonzero value for 
u k ,i, say 
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Since u kQ = we obtain u k<2 from eq 10, then u k>Zj 
etc. The fact that u ktTl+ i should turn out to be 
zero serves as a check on the calculations. 

As a very simple example, let us take the system 

u=0 at x=0, 
u=0 at x=l. 
Let h =1/6. For this system eq 6 becomes 

Vi, i+l + Vi, J-i— Vi+l, i+ Vi-1, i, 

and if we take T 7 ! =2 the array (8) turns out to 
be 

2 

10 

0-1 1 

0-101 

0—10 

0-1 

From this we have immediately 

Aq = A.2 = A± = 0, 

A 1 = A B =l i 

and the characteristic equation is 

cos 5/x+cos 3/x+cos ii = 0, 
or in factored form 

cos 3m (2 cos /x+l)=0. 
From the second form we read off the roots at once. 

_7T 2lT 37T 4w 

whence 



























1 












Ml— q> M2 — -jTf M~-Q> M4 = "^^ V5 = -^> 



57T 
6 



\ k =kTr, k = l, 2, 3, 4, 5. 

To obtain values of the characteristic functions, 
we resort to eq 10, which here becomes 

W*. 1+1=2 (C0SM*)tyfc.<— «*.«-i, 

where cos Mi=iV3, cos Ai2=i, cos ms=0, cos /z 4 = — J, 
cos /x 5 =— iV3. 



Upon calculating the w^ with this difference 
equation we have 

i = 1 2 3 4 5 6 



Wl< 


= 


] 


L V3 


2 


V3 1 


L 


^2< 


= 


] 


L 1 





-1 -] 


L 


«8* 


= 


] 


I 


— 1 


] 


L 


^4i 


= 


1 


L —1 





1 -] 


L 


%i 


= 


1 


i -Vb 


2 


-V5 ] 


I 



This particular example is remarkable in that 
the values obtained for the \ k and for the u ki not 
only are exact for the difference equation but are 
also exact for the differential system from which 
we started. We know, of course, that the correct 
X's are given by 

We readily see that the values obtained numer- 
ically for u ki are also the values of the true char- 
acteristic functions 

u x (x) = 2 sin irx, 



u 2 (x) = _ 7f sin 2tx> 
u 3 (x)=sin 3irx, 



u A (x) =— j= sin 4lttx> 

u 5 (x) =2 sin 5tx. 

When p(x) or q(x) are different from zero, the 
results will no longer be exact but only approx- 
imate. 

Several examples follow, showing the numerical 
calculations and comparing the values obtained 
with the correct values. Each of the examples 
illustrates a slightly different aspect of the process. 
Examples 1 and 2 differ only in the number of 
intervals employed. Example 3 illustrates the 
case of variable coefficients. In examples 1, 2, 3 
the coefficients satisfy the conditions of continuity 
imposed in the foregoing discussion, whereas in 
4 and 5 these conditions are not satisfied. 

The values of the characteristic functions were 
also computed in the case of Example 4, and 
compared with the correct values. 
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-> 



Characteristic numbers: Example 1 

d u dit 
DifTerential equation: -*- 2 +2 ,- -f-X 2 w = 

Boundary conditions: u = at z=0 and at x=l 
Interval: h=l/S 
Difference equation: Vij+i=jr [7Vi-i i4-9F»+i,,] — V,-,,-_i 

Boundary conditions: "^0,1 — ^8,7 = 






1 


2 


3 


4 


5 


6 


7 


8 





2097152 





























917504 























- 1064960 





802816 























-946176 





702464 

















512 





-840448 





614656 

















1120 





-746368 





537824 











748 





1764 





-662676 





470596 











1519 





2401 





-588245 









Characteristic equation: 

F 00= cos 7m+ (1.09375) cos 5/x+ (1.15869) cos 3 M + (1.19191) cos M =0 



k 


Mfc 


a M 


b A fc 


Relative error 




Degrees 








1 


23. 561 


3. 2898 


3. 2969 


-0. 00215 


2 


45. 448 


6. 3457 


6. 3623 


-0. 00261 


3 


67. 686 


9. 4507 


9. 4777 


-0.00285 


4 


90 


12. 5664 


12. 6061 


-0. 00315 


5 


112.314 


15. 6820 


15. 7398 


-0. 00367 


6 


134. 552 


18. 7870 


18. 8761 


-0. 00472 


7 


156. 439 


21. 8430 


22. 0139 


-0. 00776 



»X A . = calculated value. 
bX^ correct value. 
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Characteristic numbers: Example 2 

d il du 
Differential equation: -^ -j-2 -j — \-\*u = 

Boundary conditions: u = at x=0 and at x=l 
Interval: ^=1/16 

Difference equation: V iiJ+1 = yg [15 !%•-!,,•+ 17 V,+i,,]~ Vij-i 

Boundary conditions: V Of1 =ViQ ti =0 



^ 






1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 





115292 





















































54043 











G 



































-57871 





50665 















































-54465 





47499 









































1. 7592 





-51259 





44530 





u 



































4. 1232 





-48241 





41747 



































2. 6216 





6. 9578 





-45400 





39138 



































5. 7273 





10. 147 





-42726 





36692 





























3.4636 





9. 1926 





13.590 





-40208 





34399 





























7. 2869 





12. 910 





17. 199 





-37838 





32249 











Q 











4. 2787 





11.356 





16. 788 





20. 902 





-356C8 





30233 























8.7900 





15. 573 





20. 748 





24. 634 





-335C8 





28344 

















5.0606 





13. 431 





19. 856 





24. 723 





28. 344 





-31532 





26572 


T 














10. 225 





18. 116 





24. 136 





28. 658 





31. 985 





-29672 





24 HI 











5.8035 





15. 403 





22. 772 





28. 354 





32.508 





35. 522 





-27921 





23354 











11. 581 





20.519 





27. 339 





32. 463 





36. 233 





38. 924 





-26274 









Characteristic equation: 

F(m) = cos 15/i+ (1.05469) cos 13 M + (1.10275) cos 11 M + (1.14366) cos 9 M 

+ (1.17695) cos 7 M + (1.20225) cos 5/x+ (1.21928) cos 3m 
+ (1.22784) cos M =0. 



k 


v-u 


a x: 


b x* 


Relative error 




Degrees 








1 


11.800 


3. 2952 


3. 2969 


-0. 00052 


2 


22. 769 


6. 3583 


6. 3623 


- . 00063 


3 


33. 917 


9. 4715 


9. 4777 


- . 00065 


4 


45. 112 


12. 5976 


12. 6061 


- .00067 


5 


56. 325 


15. 7289 


15. 7398 


- . 00069 


6 


67. 546 


18. 8625 


18. 8761 


- . 00072 


7 


78. 772 


21. 9974 


22. 0139 


- . 00075 


8 


90 


25. 1327 


25. 1526 


- . 00079 


9 


101. 228 


28. 2681 


28. 2920 


- . 00084 


10 


112. 454 


31. 4030 


31. 4318 


- . 00092 


11 


123. 675 


34. 5366 


34. 5720 


- . 00102 


12 


134. 888 


37. 6679 


37. 7124 


- . 00118 


13 


146. 083 


40. 7940 


40. 8529 


- . 00144 


14 


157. 231 


43. 9072 


43. 9937 


- . 00197 


15 


168. 200 


46. 9703 


47. 1345 


- . 00348 



«> fc = calculated value. 
b X jfc = correct value. 
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Characteristic numbers: Example 3 

d 2 u 
Differential equation: 4z j~^ + X 2 u = 

Boundary conditions: u = at x=l and at x=4 

Change variable: x=s 2 

d 2 u 1 du 
Normal equation: -rj -5 — h\ 2 u = 0, u = at s=l and at s = 2 



Interval: h = l/S 



Difference equation: V it h-i*=£ [(2t + l) F,_i. ,+ (2t- l)l r ;+i. ,]- V,-. ,•- 
Boundary conditions: y 8t ,•= Vi«, ,-=0 



8 


9 


10 


11 


12 


13 


14 


15 


16 





13178880 





























6918912 























-6644352 





7233408 























-7023744 





7534800 

















10816 





-7384104 





7824600 

















20202 





-7728000 





8104050 











8263. 667 





28497 





-8057475 





8374185 











15547 





35937. 5 





-8374185 









Characteristic equation: 

F(m) = cos 7m +(1.03333) cos 5/*+ (1.05506) cos 3/*+ (1.06578) cos/i=0. 



k 


Mfc 


a x: 


h K 


Relative error 




Degrees 








1 


22. 873 


3. 1936 


3. 1966 


-0.00094 


2 


45. 164 


6. 3060 


6. 3124 


-0.00101 


3 


67. 568 


9. 4343 


9. 4445 


-0.00108 


4 


90 


12. 5664 


12. 5812 


-0. 00118 


5 


112. 432 


15. 6984 


15. 7199 


-0.00137 


6 


134. 836 


18. 8267 


18. 8595 


-0. 00174 


7 


157. 127 


21. 9391 


21. 9997 


-0. 00275 



aXfc= calculated value. 
b \ fc = correct value 
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Characteristic numbers: Example 4 
(Differential equation with singular point) 



Differential equation: t-jH — -t- + X 2 w = 

Boundary conditions: u finite at x=0, u = at x= 1 

Interval: h=l/8 



1 



Difference equation: V,-,/+i = o-- [(2i—l)Vi-i,j 

Boundary conditions: V 0l /+i = 4Vi,,— -2Vo,/— V ,,--i 









1 


2 


3 


4 


5 


6 


7 


8 


1 


























-1 


0.25 























2 


-0. 5 


0. 1875 




















-5 


1. 03125 


-0.375 


0. 15625 

















12. 125 


-2. 5625 


0. 78125 


-0. 3125 


0. 13672 














-29. 5 


6. 203125 


-1. 9375 


0. 65430 


-0. 27344 


0. 12305 











71. 6875 


-15. 09375 


4. 68896 


-1. 62109 


0. 57422 


-0.24609 


0. 11279 








-174. 25 


36. 67407 


-11.40918 


3. 92310 


-1. 421875 


0. 51782 


-0. 22558 


0. 10473 





423. 50879 


-89. 14502 


27. 72046 


-9. 54541 


3. 44104 


-1.28174 


0. 47534 


-0. 20947 






Characteristic equation: 

F(m)= : cos 8 m + 2cos 7 m+ (-0.21438) cos 6 m+ (1.57126) cos 5 M + (-0.32580) cos 4 M - 
-| (-0.35057) cos 2 //+ (1.36154^ cos ju+ (-0.17790) =0 



(1.41947) cos 3 ii 



k 


Mfc 


a x; 


b x fc 


Relative error 




Degrees 








1 


17. 212 


2. 4032 


2. 4048 


-0. 00067 


2 


39. 518 


5. 5178 


5. 5201 


-0. 00042 


3 


62. 000 


8. 6568 


8. 6537 


+ 0. 00036 


4 


84. 623 


11. 8156 


11. 7915 


0. 00204 


5 


107. 545 


15. 0161 


14. 9309 


0. 00571 


6 


131. 028 


18. 2950 


18.0711 


0. 01239 


7 


155. 614 


21. 7278 


21. 2116 


0. 02434 



* xj = calculated value. 
b X = correct value. 
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Calculation of characteristic functions for example 4 
Vk,i-i = 2jZTi f[4i cos VkWki— (2i+V)Vk,i+i) 
V k8 = 0, V k7 =l h=l/8 



k Xss \ s ^ 





1 


2 


3 


4 


5 


6 


7 


8 


1 


6. 06564 


5. 94118 


5. 54449 


4. 90869 


4. 07722 


3. 10550 


2. 05727 


1 





2 


-4. 25902 


-3. 78297 


-2. 47048 


-0. 77879 


0. 73496 


1. 61344 


1. 66115 


1 





3 


3. 78121 


2. 76528 


0. 47188 


- 1. 30444 


-1. 38768 


-0. 14450 


1. 01194 


1 





4 


-3. 81020 


- 1. 99290 


1. 01301 


1. 35254 


-0. 49927 


- 1. 13784 


0. 20836 


1 





5 


4. 27565 


1. 20726 


-1.89415 


0. 15853 


1. 27379 


- 0. 78300 


-0. 62746 


1 





6 


-5. 44631 


-0. 14517 


1. 93836 


-1. 88249 


0. 66491 


0. 71347 


- 1. 36784 


1 





7 


8. 51801 


- 1. 96624 


-0. 52633 


1. 92273 


-2. 53211 


2. 47611 


- 1. 90027 


1 






Normalized so that Vfc,o=l 



k \^ 





1 


2 


3 


4 


5 


6 


7 


8 


1 




0. 97948 


0. 91408 


0. 80926 


0. 67218 


0. 51198 


0. 33917 


0. 16486 





2 




0. 88823 


0. 58006 


0. 18286 


-0. 17257 


-0. 37883 


-0. 39003 


-0. 23480 





3 




0. 73132 


0. 12480 


-0. 34498 


-0. 36699 


-0. 03822 


0. 26762 


0. 26447 





4 




0. 52304 


-0. 26587 


-0. 35498 


0. 13104 


0. 29863 


-0. 05468 


-0. 26245 





5 




0. 28236 


-0. 44301 


0. 03708 


0. 29792 


-0. 18313 


-0. 14675 


0. 23388 





6 




0. 02665 


-0. 35590 


0. 34565 


-0. 12208 


-0. 13100 


0. 25115 


-0. 18361 





7 




-0. 23083 


-0.06179 


0. 22573 


-0.29727 


0. 29069 


-0. 22309 


0. 11740 






True values: Jo(jokX) 



^^\^ i 
k "\^ 





1 


2 


3 


4 


5 


6 


7 


8 


1 




0. 97754 


0. 91166 


0. 80679 


0. 66993 


0. 51014 


0. 33788 


0. 16421 





2 




0. 88447 


0. 57765 


0. 18368 


-0. 16840 


-0. 37274 


-0. 38424 


-0. 23137 





3 




0. 72818 


0. 13078 


-0. 33158 


-0. 35628 


-0. 03825 


0. 25859 


0. 25597 





4 




0. 52632 


-0. 24188 


-0.33776 


0. 12078 


0. 28180 


-0. 05069 


-0. 24713 





5 




0. 30138 


-0. 40077 


0. 02667 


0. 27086 


-0. 16394 


-0. 13336 


0. 21102 





6 




0. 07789 


-0. 31640 


0. 29149 


-0. 09897 


-0. 11348 


0. 21266 


-0. 15483 





7 




-0. 12066 


-0. 07480 


0. 18221 


-0. 22704 


0. 21788 


-0. 16591 


0. 08703 
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Characteristic numbers: Example 5 

Partial differential equation: - ^ = a~ ( ~ a - ) 

Boundary conditions: u finite at x = 0, u = at x=l 
Ordinary differential equation: -j—^—~ -r--{-\ 2 u = Q 



Difference equation: F»,,-+i =„• t(2i + l) f»-i./+ C2i — 1) V i+lfi ]— V iti - X 
Boundary conditions: V j= V 8ti =0 






1 


2 


3 


4 


5 


6 


7 


8 





40320 





























25200 























-27720 





29400 























-37800 





33075 

















8820 





-45937. 5 





36382. 5 

















14371. 875 





-52820 





39414. 375 











- 1634. 0625 





18604. 6875 





-59121. 5625 





42229. 6875 











-2460.9375 





22118. 90625 





-64752. 1875 









Characteristic equation: 

F(/x) = cos 7/i+ (1.64286) cos 5/x+ (1.95982) cos 3 M + (2.10044) cos m=0 



k 


Mfc 


*k 


b x fc 


Relative error 




Degrees 








1 


27. 145 


3. 7901 


3. 8318 


-0.0109 


2 


49. 044 


6. 8479 


7. 0155 


- . 0239 


3 


69. 928 


9. 7638 


10. 1735 


- . 0403 


4 


90 


12. 5664 


13. 3237 


- . 0568 


5 


110. 072 


15. 3690 


16. 4706 


- . 0669 


6 


130. 956 


18. 2849 


19. 6158 


- . 0678 


7 


152. 855 


21. 3427 


22. 7601 


- . 0623 



*\£ = calculated value. 
h \ k = correct value. 



Los Angeles, June 20, 1949. 
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